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Abstract 

A search for mixing in the neutral D meson system was performed using semileptonic 
K~e~^i' decays. The flavor of neutral D mesons at production was tagged by the charge of the 
slow pion from the decay D*'^ D^Tr~^. The measurement was performed using 140 fb~^ of data 
recorded by the Belle detector. The yield of right and wrong sign decays arising from un-mixed and 
mixed events, respectively, was determined by a fit to kinematic observables. From the number of 
signal events we derive an upper limit for the time- integrated mixing rate ro < 1.4 x 10~^ at 90% 
C.L. 

PACS numbers: 13.20.Fc, 14.40.Lb 
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I INTRODUCTION 



By considering the quark box diagrams describing — D mixing (Fig^, it is easy to 
see why the mixing rate in the charm sector is expected to be small. The amplitude for 

c 

DO 

u 

FIG. 1: Box diagram for the mixing amplitude, 
these processes can be written as 

<DVwcak|^°> oc yc:VmVcjV:^ Kimlmf) , (1) 

i,j=d,s,b 

where V^i are the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements ^ and /C is the 
loop function depending on the masses of down-like quarks. The contribution of b quarks is 
neghgible due to the small magnitude of V*fjVub- In the limit that the d and s quark masses are 
equal, the unitarity of the CKM matrix ensures a vanishing amplitude. However, the mass 
difference of d and s quarks is insufficient to result in a large mixing rate. Expectations 
for the ratio of mixed to unmixed decay rates from the box diagram are of the order of 
XQ-io _ xo~^ 0. However, due to the smallness of the short distance contributions, long 
distance effects could be important. The resulting predictions for the mixing probability 
depend on the treatment of these effects and cover a large range, reaching values as high as 

10-3 3. 

In direct searches, the signature of mixing is neutral D meson decay products with a 
charge combination opposite to that expected for a produced flavor eigenstate. Examples of 
appropriate f decay modes for these methods are K~tt'^ and K^i'^ufi^. 

If in the decay of an initially produced a K~^7c~ or K~^i~T' final state (/) is observed, this 
represents an indication of a mixing transition f. While the semileptonic decay 

modes are more difficult to reconstruct due to the presence of a neutrino in the final state, 
the hadronic modes suffer from an irreducible background of doubly Cabibbo suppressed 
(DCS) decays. 

To search for D° — mixing we used the D*^ — >• -D°7r+, ^ ^ K^e~V decay 
chain. The flavor of the produced D meson is tagged by the charge of the slow pion (vr^) 
from the D* decay. The charge of the electron determines the charm quantum number of 
the D meson at decay. If all the flnal state particles are reconstructed, signal events yield 
a sharp peak in the Am = m{7i^ K'e^u) — m{K~ v) distribution, where m{'K^ K~ e^v) 
and m{K~ v) are the invariant masses of the selected particles. In order to use the Am 
distribution to isolate the signal events, one needs to reconstruct the momentum of the 
neutrino. Since the kinematics of mixed and un-mixed decays are the same, both can be 
selected using the Am observable. 

Un-mixed decays, D*^ D^tt'^ , —>■ K~e^v, result in the "right sign" (RS) flnal 
state charge combination, Ti^K^e^ . Mixed decays, D*^ — * D^tt'^ , — > D — K^e~V, 
produce a "wrong sign" (WS) charge combination, Tr+fT+e^. Table HI lists possible charge 
combinations of flnal state particles used in the analysis. 



d,s,b 



w" 



d,s,b 



4 



Decay 


Final state 


Notation 


D*+ ^ L)07r+, D° K'e+u, un-mixed 




Right sign (RS) 


D*+^D^TT+, ^ If ^ K+e'V, mixed 




Wrong sign (WS) 


combinatorial background 




Combinatorial sign (CS) 



TABLE I: Summary of charge combinations and notations used in the analysis. 



In the following section we describe the selection of D meson semileptonic decay candi- 
dates. Section 3 describes the neutrino reconstruction and calculation of Am. A background 
description is presented in section 4 and the proper decay time distribution of mixed and 
un-mixed decays is dealt with in section 5. The succeeding section describes the fit of Am 
distributions, the obtained signal yield and the time integrated mixing rate td. In section 7 
we evaluate the systematic errors and describe several cross-checks of the measurement. 

II SELECTION OF SEMILEPTONIC D DECAYS 
Data set 

For the present analysis, data collected by the Belle detector at the center-of-mass energy 
of the T(45') resonance, corresponding to an integrated luminosity of 140 fb~^ were used. 
A detailed description of the Belle detector can be found in Simulated events were 
generated by the QQ generator and processed with a full simulation of the Belle detector, 
usmg the GEANT package 0- About 410 fb ^ of cc quark pair Monte Carlo (MC) events 
were used, together with about 164 fb~^ light quark and about 100 fb^^ BB meson pair MC 
events, which were rescaled by the appropriate factors. 

To simulate mixed D meson decays, a generic (non-mixed) MC was used with an appro- 
priate re-weighting of the proper decay time distribution. 

Selection criteria 

Electron-positron interactions with hadrons in the final state are selected from data, 
based on the neutral cluster energy, total visible energy in the center-of-mass system (cms), 
z component of total cms momentum and the position of reconstructed event vertex p, 
with an efficiency above 99%. 

The selection criteria were optimized using the simulated sample of generic u, d, s, charm 
and B meson pair events, by maximizing the efficiency for non-mixed signal events multiplied 
by the purity of the selected sample. Since the kinematic properties of mixed and non-mixed 
events are the same, such an optimization ensures optimal selection for the mixed events as 
well. 

For candidates, tracks with impact parameter with respect to the interaction point 
in the radial direction \6r\ < 1 cm and in the beam direction \6z\ < 2 cm were considered. 
This procedure suppressed badly reconstructed tracks and tracks which do not originate 
from the interaction point. A slow pion candidate was required to have a momentum of 
less than 600 MeV/c. To reduce the background from electrons misidentified as slow pion 
candidates, we required the electron identification likelihood 0] for each track, based on the 
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information from the central drift chamber, aerogel Cherenkov counter and electromagnetic 
calorimeter, to be lower than 0.1. To illustrate the 7rf selection criteria the MC efficiency 
for reconstruction of a slow pion track was found to be around 60%. Electron candidates 
were required to have a momentum p > 600 MeV/c and an electron identification likelihood 
[3| above 0.95. By that 45% of all simulated signal electrons are reconstructed. From the 
remaining tracks in an event, charged kaon candidates were chosen with momentum above 
800 MeV/c. Using the information from the time of fiight counters, central drift chamber 
and aerogel Cherenkov counters, a relative likelihood C for a given track to be a or tt"^ 
is obtained 0. Kaon candidates were selected using C^K"^) / {C^K"^) + C^n"^)) > 0.5. Using 
this selection about 50% of signal kaons were retained. 

In order to reduce the combinatorial background and the background from BB events, 
and to improve the resolution in the Am, a requirement for the cms momentum of the kaon- 
electron system, p*{Ke) > 2 GeV/c was imposed. To further suppress mesons arising 
from BB events, we required the ratio of the second and zeroth Fox- Wolfram moments R2 
to be greater than 0.2 j8|. 

Background from — > K~tt~^, with 7r+ misidentified as (RS sample) or with K~ 
misidentified as e~ and vr"*" as (WS sample), was suppressed using a requirement on the 
invariant mass of the K — e system, M{Ke) < 1.82 GeV/c^. Another source of background 
contributing to the RS and WS samples is the K^K^ decay with either of the 

charged kaons misidentified as an electron. This background is effectively suppressed by the 
requirement \M{KK) — m£,o\ > 10 MeV/c^, where M{KK) is the invariant mass of the 
kaon and electron candidate, calculated using the kaon mass for both tracks. 

An important source of background are the electrons from photon conversions, which 
could be selected as electron or slow pion candidates. Combinations where both 7r+ and e~ 
candidates are from conversions, are almost completely suppressed by requiring M{e~^e~) < 
150 MeV/c^, where the electron mass is assigned to the vr^ candidate. A further search for 
an additional of opposite charge as the semileptonic electron candidate or as the slow pion 
candidate was performed in each event. Again, events with M{e~^e~) below 150 MeV/c^ 
were rejected. For the events satisfying all the above criteria, the neutrino reconstruction 
was performed. 



Ill NEUTRINO RECONSTRUCTION 

Four-momentum conservation in an e^e~ collision implies 

Perns P-KsKe -^rest (2) 

for the signal decay, where cms stands for the center-of-mass four-momentum of the e^e~ 
system, and the index rest indicates the four-momentum of all detected particles except 
the slow pion, charged kaon and the electron candidate^^. The variable Prest is calculated 
using all the remaining charged tracks in the event with \Sr\ < 1 cm and \6z\ < 2 cm 
and all photons with energies above 100 MeV. A first approximation for the neutrino four- 
momentum Pi, is obtained using Eq. |21 The resulting Am = M{7^^K~e'^i') — M{K~e^v) 
distribution for signal events is shown in Fig. |2K as the dashed histogram. It has a peak 
at 0.145 GeV/c^, the D*^ — mass difference, and with a full width at half maximum 
(FWHM) of 55 MeV/c^. 

Two kinematic constraints were used to improve the resolution on the neutrino mo- 
mentum. First, the square of the invariant mass of selected particles was calculated using 
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Am [GeV] ^ [QeV] 

FIG. 2: a) Am distribution for simulated signal events using the missing momentum for neutrino 
reconstruction (dashed histogram) or the missing momentum together with the constraints on 
the missing mass and the mass of the D*^ (solid histogram) as described in the text, b) Am 
distribution for simulated signal events in different bins of the momentum of the kaon-electron 
system: p*{Ke) < 2.0 GeV/c (dashed histogram), 2.0 GeV/c < p*{Ke) < 3.2 GeV/c (dotted 
histogram), 3.2 GeV/c < p*{Ke) (solid histogram). The histograms are normalized to the same 
area. 

M'^i^TisKei') = (Pj, + Pn^KeY- Foi signal events the invariant mass should equal the mass 
of the D*^ meson, m£,*±. Only events with —4 GeV/c^ < MijTsKevY < 36 GeV/c^ were 
retained. For the selected events, the magnitude of Prest was rescaled by a factor x requiring 

M{7lsKeuY = (Perns - X ■ PrestY = "^D*± , (3) 

with m£,*± fixed to the nominal value of 2.010 GeV/c^ [9]. 

As a second kinematic constraint, the square of the missing mass, Pj, is used. For events 
satisfying —2 GeV/c^ < Pj < 0.5 GeV/c^, the angle a between the direction of prest and 
the direction of the TCsKe system momentum was corrected in order to yield 

ml = (Perns - E^sKe " -^rcst)^ - pl^Ke " Prcst " '^Pir^KePrest COS = 0. (4) 

The angle a was corrected by rotating jVesi in the plane determined by the vectors Prest and 

P-KsKe- 

The Am distribution obtained using the neutrino four-momentum after the use of kine- 
matic constraints is shown in Fig. as the solid histogram. The resolution is significantly 
improved, with the FWHM being about 7 MeV/c^. It also improves at higher values of 
p*{Ke), as illustrated in Fig. Eb. Using the simulated background events it has been verified 
that the described neutrino reconstruction does not induce any peaking in the background 
Am distribution. 

The MC prediction for the reconstructed Am distribution is shown in FigOl the fraction 
of signal events in the Am < 0.18 GeV/c^ region is fs = 0.620 ± 0.001. The background is 
composed of (75.6 ± 0.2)% cc events, (5.4 ± 0.1)% light quark events and (18.9± 0.1)% BB 
events. There is a small fraction (about 1%) of signal events arising from BB events. At 
this stage the efficiency for reconstructing the signal with Am < 0.18 GeV/c^ is found to be 
(5.4 ±0.1)%. 
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FIG. 3: Simulated distribution for Am with neutrino momentum reconstruction using kinematic 
constraints. The open histogram shows the contribution of the signal and the hatched histograms 
show individual contributions to the background. 




IV BACKGROUND EVALUATION AND Am DISTRIBUTION 
Background in the RS sample 

To check the MC and to avoid systematic errors arising from a possible discrepancy with 
the data, the majority of the background was described using the data. The combinato- 
rial, charge un-correlated background in the RS Am distribution can be described using 
"combinatorial sign" (CS) track combinations (Table H}. These are composed of a slow pion 
candidate, tt^, and a kaon and electron of the same charge, K'^e^. Since such combinations 
cannot arise from semileptonic decays of D mesons, they describe random track combina- 
tions giving rise to a background with a Am value less than 0.18 GeV/c^. Fig. shows a 
comparison of Am for all MC RS background events and CS data events; the former shows 
an excess of events in the region of low Am. This is a consequence of the charge-correlated 
track combinations apart from the signal itself, present in the RS sample. 

There are additional semileptonic decays with the same charge combination of final state 
charged tracks as in the — > K~e^u decay. The most important one is K'e^vK^, 
which according to the simulation represents 3| 8.8% of all reconstructed events with Am < 



0.18 GeV/c^. Other modes, like Cabibbo suppressed tt e^u, or DO ^ n^e+uK de- 

cays, with a charged pion misidentified as a kaon[16], represent an additional 1% of simulated 
events in Fig. E^. The charge correlation is also preserved if the slow pion decays into a 
muon, and the latter is taken as a slow pion candidate (0.7% of all reconstructed events). 

These decay channels have two important properties: they tend to peak at low values of 
Am and secondly, the charge of the kaon or electron candidate, together with the charge of 
the slow pion from a D*~^ decay, carries the same information on possible mixing as in the 
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FIG. 4: a) Am distribution for all simulated RS background events (dashed histogram) and for OS 
data events (points with error bars). Histograms are normalized to the same number of entries, 
b) The same comparison after associated signal and — > K~it~^tt^ decays are removed from the 
simulated RS distribution, c) The difference between the normalized Am distributions for CS data 
and the remaining simulated RS background. 



case of the selected signal decay. In the following, these decays are treated as a part of the 
signal and described as associated signal. 

Apart from the processes included in the signal, there are a number of hadronic decays 
and decays where particle misidentification might appear, contributing to a charge-correlated 
background. Using measured branching ratios for the individual processes |^, measured and 
simulated average misidentification rates and estimates of corresponding CKM elements for 
the DCS decays, we find that the most important contribution to the charge-correlated 
background in the RS sample is the decay — *■ K~7r~^7r^. This is confirmed using the MC. 
Contributions from other decays are expected to be almost an order of magnitude smaller. 
The Am distribution and the amount of the charge correlated background were evaluated 
using the MC sample. The hadronic background accounts for 0.44 ±0.01% of all events with 
Am < 0.18 GeV/c^ and is included in the CS Am distribution when extracting the yield of 
RS events. 

After removing the associated signal processes and the contribution from K^-k^-k^ 
decays from the simulated RS background, the agreement of the Am distribution with the 
CS data is good. The comparison is shown in Fig. HId and c. The x^; calculated from the 
difference of the two, has a value of 36.5 for 40 degrees of freedom. 



Background in the WS sample 

While the CS events represent a good description of the combinatorial background for the 
RS sample, the shape of the WS sample combinatorial background is different (Fig.|St')- This 
is due to different fractions of events, where the slow pion candidate is a true irf from a D*^ 
decay and either the kaon or electron candidate (or both, at least one being misidentified) 
are from the corresponding decay. The amount of this background, which will be referred 
to as the angularly correlated background, is much smaller in the WS sample. The rest of 
the background is angularly uncorrelated and represents the majority of the background 
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in the WS sample. Therefore WS background events were described by combining a 
candidate from one event with a slow pion candidate from another event {embedded pions). 
The shape of the Am distribution using the embedded pions agrees well with the MC 
predicted angularly uncorrelated WS background distribution. The comparison is shown in 
Fig. EJd andEt. The of the difference was found to be 33.1 for 40 d.o.f. An additional 
(2.04 ± 0.05)% of angularly correlated background is not represented by using embedded 
nf. The MC simulation shows that the largest contribution to this background comes 
from —>■ K'n^n^, as already established for the RS sample. The Am distribution for 
these events was obtained from simulation and added to the background described by the 
embedded slow pions. 
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FIG. 5: a) Simulated Am distribution for WS (solid histogram) and RS (dashed histogram) back- 
grounds. Both are normalized to the same area, b) Comparison of simulated WS background 
without the contribution of angularly correlated background (solid histogram) with the background 
obtained using the embedded pions (points with error bars), c) The difference of the two (normal- 
ized) distributions shown in the previous plot. 

A preliminary binned fit to the Am distribution of RS data events was performed, 
before any proper decay time selection. For the signal and associated signal the simulated 
shape of Am, Ps{Am), was assumed. Background in this sample was described using the 
CS data events with an additional contribution from — > K~7r~^7r^ events, P/^(Am). 
The fraction of the latter was fixed in the fit to the value predicted by the MC. The data 
distribution was fitted with 

VR{Am) = ATMfPsiAm) + (1 - /f )P,^(Am)) , (5) 

where the fraction of signal events (including the associated signal) and an overall nor- 
malization Afji were the free parameters. The result of the fit was = 0.732 ± 0.003. 

In order to enhance the sensitivity to mixed decays a selection based on the proper 
decay time was used, as explained in the next section. 



V PROPER DECAY TIME DISTRIBUTIONS 

It should be noted that in the present analysis no attempt is made to obtain results on the 
mixing rate from the fit to proper decay time distributions of mesons. The decay time 
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is used only to select possible mixed decays with a higher purity. To measure the mixing 
rate from the ratio of WS and RS events, one needs to evaluate the ratio of acceptances 
for mixed and non-mixed decays for decay times larger than a chosen value. The purpose 
of this section is to describe the estimation of the acceptance ratio, in which most of the 
systematic uncertainties cancel. 

The mixed events are expected to have the following proper decay time distribution ■ 

Pws(t') oc t" e-*'/^z,o , (6) 

with T£)0 = (411.7 ± 2.7) x 10~^^ s denoting the lifetime of neutral D mesons. The 
sensitivity to mixed events can be enhanced by selecting decays with a longer proper 
decay time, as the background events tend to be concentrated at lower decay times. The 
decay time t' is evaluated using the decay length from the e~^e~ interaction point to the 
reconstructed decay vertex, and the measured momentum of the meson; due to the 
shape of the KEK-B accelerator interaction region, which is narrowest in the vertical (y) 
direction, the dimensionless proper decay time t is calculated as 

Too Py CTdo 

Py is the y component of the candidate's momentum and ?/vtx and yip are the y coordinates 
of the reconstructed K — e vertex and of the interaction point, respectively. The observed 
t distribution is smeared due to the experimental resolution. Fig. ^ shows the decay time 
distribution of RS data events with Am < 0.18 GeV/c^. 




t [Xnl t [Tj t [T„] 



FIG. 6: a) Measured proper decay time distribution for the RS events (points with error bars). 
The solid line shows the results of the fit described in the text, b) Distribution of WS background 
events, obtained by combining a candidate with an embedded vr^. The solid line is again a the 
result of the fit described in the text, c) Expected decay time distribution of — > K'e^i^ decays 
(solid line) and — > K^e~T' (dashed line) calculated using the parameters obtained from 

the fit to the RS event time distribution. 

We fit the RS decay time distribution with a signal and a background term. The former is 
an exponential convolved with a sum of three Gaussians, describing the detector resolution. 
The background contribution is composed of a lifetime (exponential) and a prompt (delta 
function) component, convolved with the resolution function. The resulting value of the 
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mean t for the signal part is 1.093 ±0.044. A similar fit is performed on the WS background 
events (Fig.lHb); using the same ansatz as for the RS background. From the fitted parameters 
we deduce the expected decay time distribution of un-mixed and mixed decays, shown in 
Fig. int. To improve the sensitivity to mixed decays, we compare the expected decay time 
distribution of mixed decays with that for background (FigEb). The optimum for the 
selection criterion is found to be t > 1.5. This requirement is imposed on both RS and 
WS events. By doing that, the final ratio of WS and RS signal events depends only on the 
ratio of decay time selection efficiencies for the two classes of events, Taking the 

ratio of efficiencies instead of applying the described selection to WS events only reduces the 
systematic error arising from an imperfect description of the decay time distributions. For 
the ratio of efficiencies for t > 1.5 we obtain gunmix^gmix _ q 420 ± 0.010, where the error is 
obtained by varying all the parameters of the fit to decay time distributions by one standard 
deviation. 

Table IH] shows the composition of the simulated sample after the proper decay time 
selection. 



Component 


Fraction [%] 


Component 


Fraction [%] 


RS 




WS 




signal 


57.7 ±0.1 


cc bckg. 


73.6 ±0.1 


assoc. signal 


12.4 ±0.1 


angularly corr. bckg. 


2.18 ±0.08 


cc bckg. 


21.9 ±0.1 


uds bckg. 


12.9 ±0.1 




0.55 ±0.02 


BB bckg. 


13.5 ±0.1 


uds bckg. 


2.63 ±0.06 






BB bckg. 


4.9 ±0.1 







TABLE II: Composition of the simulated selected samples. The indented contributions are sub- 
samples of the above component. The quoted errors on the fractions are due to the MC statistics 
only. 



VI FIT OF Am DISTRIBUTION 
Signal yield 

The Am distributions for RS and WS events, after all selection requirements, including 
t > 1.5, are shown in Fig. [7| 

The data distributions are fitted using a binned fif • The RS distribution is described 
as explained in section 4. In the fit we minimize the value of 

Ni is the number of entries in the i-th bin of the RS data distribution and VR^Anii) is given 
by the number of entries of MC signal, CS and simulated i^~7r+7r° events in this bin. 

Both errors, on and on 'P/j(Amj), are included in the definition. 
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0.14 0.15 0.16 0.17 0.18 0.14 0.15 0.16 0.17 0.18 

Am [GeV] Am [GeV] 

FIG. 7: a) Distribution of Am for RS events, b) Distribution of Am for WS events. Data is 
shown with error bars, the sohd histogram shows the fit, as explained in the text. The background 
contribution is shown as the dashed histogram. 

From the fitted values of (signal fraction) and A/r (overall normalization) we obtain 
the number of signal non-mixed events iY^g = 40198 ± 329. The quoted error includes the 
statistical error of the data as well as the error due to the limited statistics of the expected 
distributions (the same holds for the fit to the WS sample, described below). The of the 
fit with 39 degrees of freedom is 69.5. The goodness-of-fit depends strongly on the amount 
of the associated signal included, due to the different widths for the Am distributions of 

— »• K~e^v and associated decay modes. This contribution is treated as one of the sources 
of the systematic uncertainty. 

For the WS sample the same description of the signal (Ps(Am)) as for the RS sample 
is used, while the background {P^ {Am)) is described using the embedded slow pions with 
the addition of simulated angularly correlated background. The result of the fit with 

Vw{Am) = Afwif^PsiAm) + (1 - f^)P^{Am)) (9) 

is = (1.7 ± 6.0) X 10-3. The quality of the fit is good: = 37.9 for 39 degrees of 
freedom. The number of mixed events obtained by the fit is A\vs = 19 ± 67. 



Time integrated mixing rate 

The time integrated mixing rate is defined as 

^ j;^v{D^ ^ e+x-,t) ^ 2 ' ^ ^ 

where the time dependent probabilities for un-mixed and mixed events are taken to be of 
the form e"*/"^ and t^e"*/"^, respectively. The dimensionless mixing parameters x and y are 
equal to x = AM/T and y = Ar/2r, where AM is the mass difference of the two neutral D 
meson mass eigenstates, AF is the difference of the decay widths and F is the average decay 
width of the two. 

The parameter r^j is obtained from the fitted WS and RS yields, 

]\T _unmix 

= «Sr ■ l=^r . (11) 
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with the ratio of efficiencies for decay time selection given in Section 5. The result is 

td = (0.20 ±0.70) X 10"^ , (12) 
where the quoted uncertainty is statistical only. 

VII SYSTEMATIC UNCERTAINTIES 
Estimation of systematic errors 

The value of the time integrated mixing rate is obtained using the ratio of WS and RS 
decays. Since the kinematics of these decays is the same, the selection efficiencies, excluding 
the acceptance of the proper decay time selection, are expected to be the same for both types 
of decays and hence cancel in the final expression (jlip . Using the MC sample of events we 
find the selection efficiency ratio to be 0.97 ± 0.05, compatible with this expectation. 

The influence of the kinematic constraints in the neutrino reconstruction on the Am dis- 
tribution of background has been verified using the MC simulation and no artificial peaking 
has been observed. The fit to the RS sample Am distribution correctly reproduces the num- 
ber of signal events. For the simulated value of = 59561 we obtain N^^ = 59193 ± 498. 
Although the numbers of RS events agree within the statistical error of the fit, the remaining 
relative difference is propagated to an absolute error on ro of 0.001 x 10~^ and is added 
to the total systematic error. The reconstruction of WS events has also been verified using 
the simulated sample. To a sample of 2.93 x 10^ non-mixed events, roughly corresponding 
to an integrated luminosity of 180 fb^^, we added 6000 mixed D meson decays, obtaining 
a sample of decays with an input value of r§^^ = 2.05 x 10~^. The ffi (Fig. |H1) with 
the same method as used for the data yields N^^ = 234 ± 85 and N^s = 48459 ± 420. 
The former can be compared to the simulated value of iV^s ~ 255. The ffited numbers of 
signal events result in = (2.0 ± 0.7) x 10~^, which very well reproduces the input value. 
Here also the relative difference between the simulated and reconstructed A\vs is taken into 
account in the final systematic error. 




Am [GeV] 

FIG. 8: Distribution of Am for simulated events including mixed D mesons (points with error 
bars), the fit (solid line) and the background (dashed line). 

The error on the ratio of decay time selection efficiencies e^fi^ix/ g™!^ -g obtained by varying 
the fitted parameters of the RS and WS sample proper decay time distributions and leads 
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to an additional 0.005 x 10~^ error on td. We evaluated the same efficiency ratio using the 
MC sample and obtained _ q 4qq j- 0.012, which is consistent with the value 

obtained from the proper decay time fit. The proper decay time selection was additionally 
verified by performing the fit to RS and WS Am distributions using several different decay 
time requirements, t > 1.2 and t > 1.7. The resulting values of r£, are (0.30 ± 0.70) x 10^^ 
and (0.21 ± 0.70) x 10"'^, respectively. We take the largest difference as a systematic error 
due to an imperfect description of the decay time distribution. 

The next possible source of systematic uncertainty arises from MC samples used to fit 
the Am distributions. Using the measured branching fractions and corresponding errors |^ 
of decay channels included in the associated signal, and their relative contributions to the 
selected sample as obtained by MC, we estimate that the amount of this component in the 
fit needs to be varied by 40%. The fit was repeated and the difference in the obtained values 
of td is found to be 0.015 x 10"^. 

The fraction of the angularly correlated WS background is estimated from MC to be 
(2.04 ± 0.05)%. Since the majority of these events arise from — > K~tt'^7c^ decays, we 
repeated the WS fit by changing the amount of this background in accordance with the 
relative error on B[D^ — > K'ti^-k^) j^. The variation of iVws resulted in a 0.05 x 10~^ error 
on r£). 

The systematic uncertainties are summarized in Table IIIII By summing the individual 
contributions in quadrature, we obtain the total systematic error of 0.11 x 10~^. It is clear 
that the uncertainty of the result on the time integrated mixing rate is limited by the 
statistics of the data. 



Source 


Estimated by 


Resulting uncertainty [xlO^] 


u reconstruction 




±0.001 






±0.016 


ratio of t selection 


±1(7 of t distribution fit parameters 


±0.005 


efficiencies 


using different t selection values 


±0.100 


amount of associated signal 


±40% 


±0.015 


amount of charge correlated 


±7% 


±0.050 


WS backg. (L>o ^ K~t:+'k^) 






Total 




±0.11 



TABLE III: Summary of systematic measurement uncertainties. The evaluation of individual 
contributions are described in the text. 



Checks of the analysis method 

Apart from the checks on the quality of background description described in previous 
sections, several other tests of the analysis method were performed. 

The fits to the RS Am distribution were done individually for different data taking 
periods. The resulting signal fractions show no significant deviation from the overall result 
used in the analysis. 

In order to check for a possible bias in the result, we performed a fit to the Am distribution 
of WS events using an inverted decay time selection, i.e. fitting the events with t < 1.5. A 
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possible mis- assignment of the charge of the slow pion candidate, for example, would appear 
as a peak in this sample and could bias the fitted number of WS events. Since t < 1.5 is 
the region where the sensitivity to mixed events is reduced, we expect no such contribution 
and any positive signal would be an indication of a peaking background behavior. The fit 
was done in the same manner as the fit to t > 1.5 distribution and the resulting number of 
mixed events is A^ws = —77 ±88, in accordance with the expectation. The fitted number of 
non-mixed decays is N^^g = 82943 ± 456, which taking into account the corresponding ratio 
of decay time selection efficiencies gives = (—3.1 ± 3.6) x 10^^. 

By repeating the fit procedure without the decay time selection, we obtain A\vs = —48 ± 
111 and A^Rs = 123054 ± 562. While the central value of = (-0.39 ± 0.90) x lO'^ 
is compatible with the null result, as for the events with t > 1.5, the statistical error is 
increased, confirming the usefulness of the decay time selection. 

To test the bias as well as the obtained statistical error of A^ws we divided the available 
simulated WS sample into two parts. The corresponding Am distributions were then fitted 
one to another. The result of the fit is A^ws = —66 ± 62. The statistical error of the 
"background" distribution was not included in this fit. Afterward we independently varied 
the contents of each bin within the statistical error and repeated the fit many times. The 
distribution of resulting A^'ws is centered at the starting value (due to a statistical fluctuation 
when arbitrarily dividing the MC sample into two parts) with a width of 61.9 ± 0.5. The 
latter represents the expected statistical error of the fit on a sample with about the same 
statistics as the data itself. The expected error is in agreement with the fit to the data when 
not taking into consideration the statistical error of the expected distributions (±64). 



VII CONCLUSIONS 

We have searched for — D'^ mixing in semileptonic D meson decays. The resulting 
time integrated mixing rate is 

= (0.20 ±0.70 ±0.11) X 10-3 , (13) 

where the first error represents the statistical and the second the systematic uncertainty. 

Since no significant signal for the WS decays K^e~u is observed, we put an upper 

limit on the time integrated mixing rate. We combine statistical and systematic errors in 
quadrature and asume a Gaussian distribution of the resulting total error; using the Feldman 
Cousins approach 3] to estimate the upper limit in the vicinity of the physical boundary 



{fD > 0) we find 

td < 1.4 X 10-3 at 90% C.L. (14) 

The limit obtained is more restrictive than the existing world average 

T{K+tv)/T{K-l+v) < 5 X 10-3ja]. Using expression ^ and the value of 
Ar/r = 2y = 0.016 ± 0.010j|£], the above result can be expressed as 

AM < 13 X 10^° hs-^ at 90% C.L. (15) 
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